Macroscopic entanglement by entanglement swapping 
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We present a scheme for entangling two micromechanical oscillators. The scheme exploits the 
quantum effects of radiation pressure and it is based on a novel application of entanglement swapping, 
where standard optical measurements are used to generate purely mechanical entanglement. The 
scheme is presented by first solving the general problem of entanglement swapping between arbitrary 
bipartite Gaussian states, for which simple input-output formulas are provided. 
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Micromechanical resonators with fundamental vibra- 
tional mode frequencies in the range 10 MHz— 1 GHz can 
now be fabricated Q. Applications include fast, ultra- 
sensitive displacement detectors |2(, electrometers, and 
radio frequency signal processors. Advances in the de- 
velopment of micromechanical devices also raise the fun- 
damental question of whether mechanical systems con- 
taining macroscopic numbers of atoms will exhibit quan- 
tum behavior. In particular, it is interesting to see under 
which conditions it is feasible to prepare micromechan- 
ical oscillators in entangled states, where the quantum 
nature becomes most manifest. Proposals of this kind 
recently appeared in the literature. One could entan- 
gle a nanomechanical oscillator with a Cooper-pair box 
[3, with an ion 0, or with single photons moreover 
Ref. studied how to entangle an array of nanomechan- 
ical oscillators, Ref. [jj proposed to entangle two mirrors 
of an optical ring cavity, while Ref. Q considered two 
mirrors of two different cavities illuminated with entan- 
gled light beams. In this Letter we will show how to en- 
tangle two micromechanical oscillators, which are parts of 
two different micro-opto-mechanical systems, by means 
of continuous variable (CV) entanglement swapping per- 
formed on the "optical parts" . Our proposal exploits 
some results of Ref. 9] , which showed how radiation pres- 
sure can entangle a vibrational mode of a mirror with 
the back-scattered sideband modes of an intense optical 
beam. Thanks to this optomechanical entanglement, one 
can entangle two different mechanical oscillators by per- 
forming appropriate optical measurements. Note that, 
very recently, Ref. ^(J nas experimentally observed the 
two considered sideband modes, although no measure- 
ments of entanglement have been performed. The pro- 
posal is presented by first solving the general problem of 
CV entanglement swapping between arbitrary Gaussian 
bipartite states, for which we are able to derive compact 
input-output formulas. Such a study suitably generalizes 
the previous fully optical theoretical analyses and it 
has a straightforward application to our physical scheme. 

Let us first study entanglement swapping with Gaus- 
sian states in general. We consider a pair of CV systems 



labelled by k = 1,2, where each of them is described by 
a pair of conjugate dimensionless quadratures Xk and pk- 
Introducing the vector £ T = (ii,pi, X2,P2), we can write 
the canonical commutation relations in the compact form 
[£>lm] =iJim (l,m= 1,...,4), where 
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and © denotes the usual direct sum operation. An arbi- 
trary state pu of the two oscillators can be described by 
its Wigner characteristic function $12 (0 = Tr[pi 2 l)(C)], 
where D(Q = exp(— i£ T £) is the Weyl operator and 
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Such a state is said to be Gaussian if the 



corresponding characteristic function is Gaussian, i.e., 
$12 (C) — exp(— ^ r V^/2+ia r Q). In such a case, the state 
P12 is fully characterized by its displacement d = (£) 
and its correlation matrix (CM) V, whose generic ele- 
ment is defined as Vi m = (A£;A£ m + A£ m A£;)/2 where 
A£/ = — (£;). All the information about the quan- 
tum correlations between the two oscillators is contained 
in the CM, which is a 4 x 4, real and symmetric ma- 
trix, satisfying the uncertainty principle V + > 0. 
Putting the CM into the form 



V=\ A D 
1 1 D T C 
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where A, C, and D are 2x2 real matrices, one can easily 
derive a quantitative measure of entanglement. In fact, 
from the CM J5J, one can compute the two quantities 



+ = o-l/2 
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where S(V) = det A + det C — 2 det D. Such quantities 
correspond to the symplectic eigenvalues of the ma- 
trix AfcVAj,, which is obtained from V through the partial 
transposition (PT) transformation A& = A[l, 1,1,-1], 
where A[Ai, A„] denotes a diagonal matrix with the 
n entries Xk in the diagonal. One can then prove [l2T | 
that the minimum PT symplectic eigenvalue rf~ repre- 
sents an entanglement monotone, since it is connected to 



the logarithmic negativity by 
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Ejsf — max[0, — In 2^ ] 
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In particular, the Gaussian state p\ 2 is entangled if and 
only if 77 < 1/2, which is equivalent to > accord- 
ing to Eq. Q). 

In general, CV entanglement swapping involves two 
pairs of modes, (a, ci), owned by Alice and Charlie, and 
(6,02), owned by Bob and Charlie. We assume that the 
two pairs are described by two arbitrary Gaussian states 
p ac and pbc, having generic CMs 
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where A, B, C\, C 2 , Dt, D 2 are 2x2 real matrices. The 
basic idea of entanglement swapping is to transfer the bi- 
partite entanglement within the near pairs Alice-Charlie 
and Charlie-Bob to the distant pair Alice-Bob, by means 
of a suitable local operation and classical communica- 
tion performed by Charlie. In fact, Charlie mixes his two 
modes c\ and c 2 through a balanced beam splitter (BS) 
and then he measures the output quadratures x C2 — x ci 
and p C2 + p ci . After this measurement, which realizes 
a CV version of a Bell measurement, he communicates 
the outcomes to Alice and Bob [ll]. The final output 
state pah of Alice and Bob turns out to be a Gaussian 
state whose CM can be expressed in terms of the input 
matrices' blocks as 



Vab = Vll 
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where V a % 



A® B, A a = R® I and A b = / © R are 
the PT transformations over Alice and Bob respectively, 
with R = A[l,-1] and I = A[l,l], and finally g' 1 = 
detCi 4- detC 2 + Tt{RC 2 RJC\J t ). The general result 
of Eq. © specializes to a very simple form if the input 
matrices V a b and Vb c are put in the normal form, i.e., A — 
al, B = 67, Cfe = Cf.1 and D/, = A[dk,d' k ] for k = 1,2. 
This can always be achieved by means of local symplectic 
transformations 01 a * Alice, Bob and Charlie's sites, 
which therefore do not affect the entanglement resources 
of the initial states p ac and pb c - In such a case, we have 



V ab = 



al 



bl 
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A[d?,di 2 ] A[-did 2 ,di4] 
A[-did 2 ,did 2 ] A[dld' 2 2 } 



providing a straightforward expression for the CM of the 
final Alice and Bob's state. 

In general, the output entanglement cannot be ex- 
pressed as a simple function of the input entanglement. 




FIG. 1: Two identical intense laser pulses are incident on two 
identical micromechanical oscillators. Each beam excites a vi- 
brational mode, which generates two optical sideband modes 
in the back-scattered field. Both reflected Stokes modes are 
mixed through a balanced BS and the output ports are ho- 
modyned. The anti-Stokes modes can be either ignored (non- 
assisted strategy) or heterodyned (assisted strategy). 



Nonetheless, this becomes possible in the symmetrical 
case of identical physical resources shared by the two 
pairs, i.e., when p ac = pb c = pi n . If the input matrix 
Vi n of pi n is converted from the general form of Eq. @ 
to its normal form, the transformation connecting Vi n to 
the output matrix V ou t of p ou t = p a b can be written as 



V m = 



al A[d,d'] 
A[d,d'] ci 

( A[2ac- d 2 ,2ac- 
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Note that in Eq. iJHJ the output state is a symmetric 
Gaussian state, whose entan glem ent properties have been 
widely studied in literature 14]. By applying Eq. to 
Vi n and V ou t, one can connect the PT symplectic eigen- 
values r/f ut of the output matrix with the ones rjf^ of the 
input matrix by the simple input-output relations 



Vout = c VlnVi 

Vout = a ■ 
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Note that the output eigenvalues are uniquely determined 
by the local symplectic invariants of the original CM of 
Eq. J5J, since here a = V det A and c = v det C. 

Let us now apply the above results and, in particu- 
lar, the input-output relation of Eq. @, to see how one 
can entangle two micromechanical oscillators by means of 
entanglement swapping. Consider the system of Ref. 0, 
which is schematically depicted in the left part of Fig. ^ 
The radiation pressure of an intense monochromatic laser 
beam (frequency uj), incident on a micromechanical os- 
cillator, generates an effective coupling between a vibra- 
tional mode of the oscillator (label a, frequency 51) and 
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FIG. 2: Output log-negativity Effi versus time i, for T = OK 
(solid), T = 5mK (dotted), T = 300K (dashed) (see text for 
the other parameters), (a) refers to the non-assisted strategy; 
(b) refers to the assisted strategy. 

the two first optical sideband modes induced by the vi- 
brations in the back-scattered field, i.e., the Stokes mode 
(label c, frequency u> — ft) and the anti-Stokes mode (la- 
bel d, frequency w + fi) (e.g., see Ref. 0, where these 
sidebands have been generated and detected). The vi- 
brational mode is characterized by a relaxation time T _1 
which globally takes into account several damping effects, 
like internal losses of the medium and those due to the 
clamping. All these damping effects can be neglected 
if the duration of the laser pulse is much shorter than 
r _1 so that the system dynamics can be assumed uni- 
tary during the interaction time. In this unitary descrip- 
tion of the process, the interaction Hamiltonian for the 
vibrational mode, the Stokes mode and the anti-Stokes 
mode is given by H = —ihx(ac— a^c^) — ihO{a) c' — ac't), 
where \ an d are coupling constants whose ratio r = 

1/2 

9/x= [(<a> + fi)/(w — Q)] > 1 depends only on the in- 
volved frequencies, and the operators d, c and & are the 
annihilation operators of the three modes a, c and d, 
satisfying [o, aH = \c,c'] = 1, etcetera. Starting from 
vacuum states for the optical modes and a thermal state 
for the vibrational mode, with mean excitation number 
n = [exp(Ml/KT) — (T being the temperature), the 
evolved state is Gaussian with the 6x6 CM 

/ (A +1/2)1 CR TR \ 

Vcad = CR (B + 1/2)1 -VI 

\ TR -VI (£ + 1/2)1 J 

(11) 

where A, B, C, V, £ and T, explicitly given in Ref. 
depend on r, n and the scaled time t' = t(9 2 — X 2 ) 1 ^ 2 - 

Suppose that Alice and Bob are supplied with two 
identical micromechanical oscillators, equally illuminated 
by an intense laser pulse, and suppose that the reflected 
Stokes and anti-Stokes modes are sent to Charlie (i.e., to 
our experimental detection apparatus). In a first strat- 
egy (see Fig.^l, Charlie simply traces out the anti-Stokes 
modes and performs a CV Bell measurement over the re- 
maining Stokes modes. In such a case, just after reading 



out the measurement results at time t, the two mechan- 
ical oscillators are described by a Gaussian state p out 
with CM V out given in Eq. |JS} through the replacements: 
a =B + 1/2, c = A+ 1/2 and d = -d' = C. By applying 
Eq. one can compute the entanglement r]~ ut of the 
output state. Adopting the parameters (# 2 — x 2 ) 1 ^ 2 — 10 3 
s _1 , r — 1 + 2.5 x 10~ 7 and fl = 5 x 10 8 s _1 , we have stud- 
ied the output log-negativity Eff* — max[0, — ln2?7~ ut ] 
as a function of the interaction time t, and for differ- 
ent values of the external temperature T. In Fig. Efa) 
one can see the existence of wide non-zero time regions 
giving a strong entanglement between the two mechan- 
ical modes, even for large temperatures. For instance, 
if the Bell measurement is performed after an interac- 
tion time t ~ 0.8/is, one gets an optimal entanglement 
measure Effi ~ 0.88, which does not depend on tem- 
perature (see the maximum in Fig. |2fa)). This is an 
effect of the quantum interference which holds provided 
that the mechanical damping time T _1 is much longer 
than the duration of the laser pulse and therefore the 
unitary treatment of the process is valid. However, the 
living time of the generated entanglement will depend 
on the damping time T _1 scaled by a factor depend- 
ing on temperature. This decoherence time is given by 
7- 1 = T- 1 ln[(2n + 1 - 2 Vmt )(2n + e" 1 - 2e~Sou t)~\ 
which gives the limits 7 1 ~ 1/V for T — ► and 
7- 1 ~ hn{2TKT)- 1 for KT > Ml. Therefore, even for 
long r _1 , the entanglement living time will be affected 
by temperature and it will be very short for high temper- 
atures. For instance, if Q, ~ 5 x 10 s s" 1 and ~ Is, 
one has 7 -1 ~ 3 fis at T ~ 300 K, which has to be com- 
pared with 7" 1 ~ 1 ms at T ~ 1 K. It is therefore nec- 
essary to consider low temperatures in order to preserve 
the mechanical entanglement, even if this assumption is 
not crucial for its generation. Now, since the effective 
masses of the vibrational modes can be of the order of 
10~ 10 kg, this result proves the feasibility of entangling 
massive macroscopic objects combining radiation pres- 
sure effects and simple optical homodyne measurements. 
Actually, this performance can be further improved by 
assisting the Bell measurement on the Stokes modes by 
additional heterodyne measurements on the anti-Stokes 
modes (see Fig. [Q. In fact, heterodyne measurements 
are the Gaussian conditional measurements optimizing 
the en tang lement between the oscillator and the Stokes 
mode 15]. Adopting this modified protocol, the out- 
put state p ou t will be again Gaussian, but with a modi- 
fied CM V ou t, as given by Eq. JHJ with the replacements: 
a = B + 1/2 - V 2 /{£ + 1), c = A + 1/2 - T 2 /{£ + 1) 
and d = —d! = C + VT /(£ + 1). The improvement is 
shown in Fig. EJb) where the entanglement time regions 
are wider and the maximum achievable entanglement is 
larger than before (Effi ~ 1.1 at t ~ 1/xs). 

An important related issue is the experimental de- 
tection of the generated macroscopic entanglement. In 
general, proving the existence of entanglement between 
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two oscillators requires the measurement of two differ- 
ent linear combinations of quadratures, like the rela- 
tive distance X r< a = x a — Xf, and the total momentum 
P tot = p a + p b 0, 0. However, in the present case, 
the final Gaussian state of Alice and Bob, both in the 
assisted and non-assisted case, is a symmetrical channel 
where the variance of relative distance and total momen- 
tum are equal, i.e., ((Ai rel ) 2 ) = ((AP tot ) 2 ) (from Eq. © 
with d = —d'), and, in particular, ((Al re i) 2 ) = 2r}~ ut 
(from Eq. ©). If the two micromechanical oscillators are 
highly reflecting mirrors, they can be used as end mirrors 
of a Fabry-Perot cavity. It is known that when this cavity 
is resonantly driven by an intense laser field, the detec- 
tion of the phase quadrature of the cavity output pro- 
vides a real-time quantum non-demolition measurement 
of the mirror relative distance X re i 16]. Thus, if the 
cavity is driven soon after Charlie's measurements and 
one performs a homodyne detection of the output field, 
one can measure the variance ((AX re i) 2 ) and, therefore, 
the log-negativity of the swapped state via Eq. Q. If 
the two oscillators are different, the final state of Alice 
and Bob is no more symmetrical and one has to mea- 
sure both ((AX re i) 2 ) and ((AP (of ) 2 )in order to surely 
detect entanglement (e.g., see Ref. |l7|). Alternatively 
one can measure x and p of each oscillator by shining 
a second, intense "reading" laser pulse on it, and ex- 
ploiting again the same optomechanical scattering pro- 
cess. In fact, as shown in Ref. by performing a het- 
erodyne measurement of the linear combination of the 



two sidebands 



d\ one can reconstruct the oscillators' 



state. Therefore one can detect the generated mechan- 
ical entanglement provided that the living time of the 
entanglement, 7 -1 , is long enough and this is achiev- 
able at cryogenic temperatures. In fact, as we have seen 
above, the thermal environment makes the mechanical 
entanglement decay with a relaxation time which can be 
quite long at low enough T (7 ~ 1 ms at T ~ 1 K if 



n ~ 5 x 10 8 s- 1 and r- 



ls). 



swapping. The scheme exploits the optomechanical en- 
tanglement between each oscillator and the reflected side- 
band modes of an intense laser field, which can be gen- 
erated by radiation pressure. Optical measurements are 
then able to swap this entanglement to the oscillators, 
with the non-trivial effect of changing the entanglement 
from optomechanical to purely mechanical. To study the 
scheme, we have first solved the general problem of en- 
tanglement swapping between generic bipartite Gaussian 
states, providing simple input-output formulas. 
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In conclusion, we have presented a scheme for entan- 
gling two micromechanical oscillators by entanglement 
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